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A theoretical attempt to identify the physical process responsible for the afterglow emis- 
sion of Gamma-Ray Bursts (GRBs) is presented, leading to the occurrence of thermal 
emission in the comoving frame of the shock wave giving rise to the bursts. The deter- 
mination of the luminosities and spectra involves integration over an infinite number 
of Planckian spectra, weighted by appropriate relativistic transformations, each one 
corresponding to a different viewing angle in the past light cone of the observer. The 
relativistic transformations have been computed using the equations of motion of GRBs 
within our theory, giving special attention to the determination of the equitemporal 
surfaces. The only free parameter of the present theory is the "effective emitting area" 
in the shock wave front. A self consistent model for the observed hard-to-soft transition 
in GRBs is also presented. When applied to GRB 991216 a precise fit (x^ ~ l.OTs) of 
the observed luminosity in the 2-10 keV band is obtained. Similarly, detailed estimates 
of the observed luminosity in the 50-300 keV and in the 10—50 keV bands are obtained. 

Keywords: black hole physics — gamma rays: bursts — gamma rays: observations — 
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1. Introduction 

Gamma- Ray Bursts (GRBs), following the observations by the BATSE instrmnent 
(Paciesas et al.Ql), have been characterized by a few global parameters (see e.g. 
Fishman & Meega such as the fluence, the characteristic duration (Tgo or T50), 
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and the global spectral distribution given e.g. by the Band relation (Band et al\^. 
After the discovery of the afterglow (Costa et al.l^, two additional important pa- 
rameters have been added: the power-law indexes of the afterglow and the source 
luminosity. 

It has become clear that a variety of different eras are present in the GRB 
data and that GRBs quite possibly have the most extreme time-variation of any 
phenomena in nature (see e.g. RufRni et al0 and references therein). We present 
here an attempt to derive from first principles the instantaneous GRB luminosity in 
selected energy bands and the GRB spectra. We use GRB 991216 as the prototype 
(see RufRni et a\^ since this source offers a superb set of data by BATSE in the 50- 
300 keV band (see BATSE Rapid Burst ResponsJH) and by R-XTE and Chandra 
in the 2-10 keV band (see Piro et al.EI, Corbet & SmitlP) to be compared to the 
theoretically predicted ones in the 2-300 keV range. We also give physical reasons 
for the often mentioned hard-to-soft transition observed in the majority of GRBs 
(see e.g. Frontera et al.^l Ghirlanda et al.HI, PirarU^, Piro et al!^). 



2. The model 

The complete dynamics of GRB 991216 has been computed (see Ruffini et alP|). 
The initial conditions we adopted for this source at t = 10~^^ s s are a spherical 
shell of e"'"-e~-photon neutral plasma lying between the radii tq = 6.03 x 10^ cm and 
ri = 2.35 X 10^ cm: the temperature of such a plasma is 2.2 MeV, the total energy 
Etot = 4.83 X 10^3 erg and the total number of pairs A^e+e- = 1-99 x 10^^. These 
conditions have been derived from vacuum polarization processe s oc curring in the 
dyadosphere of an ElectroMagnetic Black Hole (EMBH) (RuffinPI, P reparata et 
aM). 

The optically thick electron-positro n pl asma self-propels itself outward reaching 
ultrarelativistic velocities (Ruffini et alM^, then interacts with the remnant of the 
progenitor star and by further expansion becomes optically thin (Ruffini et al.*^). 
As the transparency condition is reached, the Proper-GRB (P-GRB) is emitted with 
an extremely relativistic shell of Accelerated Baryonic Matter (the ABM pulse, see 
Ruffini et al.CSl). 

It is this ABM pulse which gives origin to the afterglow through 
its interaction with the ISM, whose average density is assumed to be (riism) = 1 
particle/cm'^. In such a collision the "fully radiative condition" is implemented (see 
Ruffini et al.'^for details): the internal energy Ai?int which results is instantaneously 
radiated away. 

The equations of motion in our model depend only on two free parameters: 
the total energy Etot, which coincides with the dyadosphere energy Edya, and the 
amount Mb of baryonic matter left over from the gravitational collapse of the pro- 
genitor star, which is determined by the dimensionless parameter B — Mbc^ / E^ya- 
These two parameters have been determined by fitting, with high accuracy, the 
bolometric intensity and the slope of the afterglow (Ruffini et aL^i^ . We have also 
fit, again with high accuracy, the substructures observed in the E-APE which result 
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from inhomogeneities in the ISM, stil l ma intaining; an average density distribution 
(nism) — 1 particle/cm^ (Ruffini et alP^. 

3. The newly assumed origin of the afterglow X- and 7-ray 
radiation 

Here we adopt three basic assumptions: a) the resulting radiation as viewed in 
the comoving frame during the afterglow phase has a thermal spectrum and b) 
the ISM swept up by the front of the shock wave, with a Lorentz gamma factor 
between 300 and 2, is responsible for this thermal emission. We also adopt, like 
in our previous papers (Ruffini et al.'^'-' '-' ^^), that c) the expansion occurs with 
spherical symmetry. These three assumptions are different from the ones adopted 
in the GRB literature, in which the afterglow emission is believed to originate from 
synchrotron emission in the production of the shock or reverse shock generated 
when the assumed jet-like ejecta encounter the external medium (see e.g. Giblin et 
al.l^ and references therein). 

The structure of the shock is determined by mass, momentum and energy con- 
servation, i.e., the constancy of the specific ent halpy, which are standard conditions 
in shock rest frames (Zel'dovich & Rayze J22I) and have already been used in our 
derivation (Ruffini et al.l^. The only additional free parameter of our model is the 
size of the "effective emitting area" in the shock wave front: A^ff. 

The temperature T of the black body in the co-moving frame is then 

T=( ^^'"^ V^" (1) 

where 

7^ = 4^ (2) 

is the ratio between the "effective emitting area" and the ABM pulse surface Aabm, 
AE'int is the internal energy developed in the collision with the ISM in a time 
interval At in the co- moving frame and a is the Stefan-Boltzmann constant. The 
ratio TZ, which is a priori a function that varies as the system evolves, is evaluated 
at every given value of the laboratory time t. 

All the subsequent steps are now uniquely determined by the equations of mo- 
tion of the system. The basic tool in this calculation involves the definition of the 
EQuiTemporal Surfaces (EQTS) for the relativistic expanding ABM pulse as seen 
by an asymptotic observer. See Fig. 1 in Ruffini et aL^^and Bianco & Ruffin^^^ 2^ 

We are now ready to evaluate the source luminosity in a given energy band. The 
source luminosity at a detector arrival time f^, per unit solid angle dil and in the 
energy band [vi, 1^2] is given by (see Ruffini et al.El): 

dEl;"''"'^ f Ae „ dt 



dtidn 



f ^ „ cos i9A-4 -^1^(^/1,1.2, T,„)dS, (3) 

J EQTS ^'^ "''o 



4 R. Ruffini, C.L. Bianco, P. Chardonnet, F. Fraschetti, V. Gurzadyan, S.-S. Xue 



where Ae — AEint/V is the energy density released in the interaction of the ABM 
pulse with the ISM inhoniogeneities measured in the comoving frame, A = 7(1 — 
{v I c) cos ■&) is the Doppler factor, W (yx , , Tarr) is an "effective weight" required to 
evaluate only the contributions in the energy band [i^i , 1^2] , is the surface element 
of the EQTS at detector arrival time on which the integration is performed (see 
also Ruffini et alPJ) and Tarr is the observed temperature of the radiation emitted 
from dS: 

_ T 1 

- 7(1- icosd) (TT^ ■ 

The "effective weight" W [y\^ Tarr) is given by the ratio of the integral over 
the given energy band of a Planckian distribution at a temperature Tarr to the 
total integral aTa„: 



W (vi,V2,Tarr) = ^f;^ I p(Tarr,v)d\ 1 , (5) 

O-^ arr J ui \ ^ / 

where p {Tarr, i^) is the Planckian distribution at temperature Tarr- 

2 hv 

= (6) 

4. The best fit of the observed flux in selected energy bands 

We can now proceed to the best fit of the observed data using GRB 991216 as 
the prototype. Such an estimate is perfectly well defined from a theoretical point 
of view, although from a numerical point of view the integration on all the EQTS 
and the associated relativistic transformations have raised unprecedented and time 
consuming difficulties. Almost 10^ paths with different temperatures and dilferent 
Lorentz boosts had to be considered in the integration over the EQTS. We give 
in Figs. EHl the results for the three energy bands 50-300 keV (BATSE), 2-10 
keV (R-XTE, Chandra) and 10-50 keV. It is most remarkable that the best fit is 
obtained simply by a factor TZ, which is monotonically varying in the range: 

3.01 X IQ-^ > 7^ > 5.01 X 10^12 (7) 

respectively in correspondence with the beginning of the afterglow emission and 
the last observation by Chandra at ~ 37 hr after the GRB. We point out the 
perfect agreement with the data obtained by BATSE Rapid Burst ResponsJ^ in 
the energy range 50-300 keV (see dashed line in Fig.^). Very impressive is the fit 
of the data obtained by the R-XTE and Chandra satellites (see Halpern et al.E^) in 
the energy range 2-10 keV (see dotted line in Figs.n]-|2Jl. These data are fitted with 
a ~ 1.078. This fit covers a time span of ^ 10^ s and is impressive if we recall 
that it is a function of the single parameter TZ. The fit can be further improved, 
reaching a ~ 0.48, when a radial dependence in (riism) is introduced, ranging 
from (riism) — 1 particle/cm"^ in the E-APE region (r ~ 5 x 10^^ cm) to (riism) ~ 3 
particle/cm"^ in the latest afterglow phases (r ~ 4 x 10^^ cm). Both in Fig. ^ and 
Fig. 13 the solid fine gives the bolometric luminosity (see details in Ruffini et al."^. 
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Fig. 1. Best fit of the afterglow data of GRB 991216. The sofid curve is the bolometric luminosity. 
See RufHni et al.^** for the radial approximation and RufBni et al.^'^^ for the relativistic analysis 
of the off-axis contributions. The three dotted curves correspond to the luminosities in the bands 
50-300 keV, 10-50 keV and 2-10 keV respectively. Near the E-APE, where the BATSE data are 
present, almost all the luminosity is in the 50-300 keV band. The afterglow data from R-XTE and 
Chandra (see Halpern et al.^^) in the 2—10 keV are perfectly fit by the corresponding luminosity 
curve (see also Fig. 1^. 

5. On the time integrated spectra and the hard-to-soft spectral 
transition 

We turn now to the much debated issue of the origin of the observed hard-to- 
soft spectral transition during the GRB observations (see e.g. Frontera et al.^^, 
Ghirlanda et al.^^ PirarU^, Piro et al.U^ . We consider the instantaneous spectral 
distribution of the observed radiation for three different EQTS: 

• = 10 s, in the early radiation phase near the peak of the luminosity, 

• — 1.45 X 10^ s, in the last observation of the afterglow by the Chandra 
satellite, and 

• = 10"* s, chosen in between the other two (see Fig.O. 

The observed hard-to-soft spectral transition is then explained and traced back to: 



(1) a time decreasing temperature of the thermal spectrum measured in the co- 
moving frame, 
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Fig. 2. This is an enlargement of Fig. 0in tlie region of tiic afterglow data in the 2-10 kcV band 
from the R-XTE and Chandra satellites, showing the perfect agreement between the theoretical 
curve and the observational data. The reduced value for this fit is — 1.078. 

(2) the GRB equations of motion, 

(3) the corresponding infinite set of relativistic transformations. 

A clear signature of our model is the existence of a common low-energy behavior of 
the instantaneous spectrum represented by a power-law with index a = +0.9. This 
prediction will be possibly verified in future observations. 

Starting from these instantaneous values, we integrate the spectra in arrival time 
obtaining what is usually fit in the literature by the "Band relation" (Band et alP^. 
Indeed we find for our integrated spectra a low energy spectral index a — —1.05 
and an high energy spectral index /3 < — 16 when interpreted within the framework 
of a Band relation (see Fig.^]). This theoretical result can be submitted to a direct 
confrontation with the observations of GRB 991216 and, most important, the entire 
theoretical framework which we have developed can now be applied to any GRB 
source. The so obtained theoretical predictions on the luminosity in fixed energy 
bands can be then straightforwardly confronted with the observational data. 

6. Conclusions 

In addition to the above results, we have also applied our model to GRB 980425, 
which is one of the weakest GRBs observed, with an energy of the order of ^ 10**^ 
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Fig. 3. The instantaneous spectra of the radiation observed in GRB 991216 at three different 
EQTS respectively, from top to bottom, for = 10 s, = 10* s and t^ = 1.45 X lO'' s. These 
diagrams have been computed assuming a constant (nism) — 1 particle/cm'^ and clearly explains 
the often quoted hard-to-soft spectral evolution in GRBs. 



ergs (see Ruffini et al i'^^ l ^^l RufSnPS). Our model then applies over a range of 
energies spanning 6 orders of magnitude. 

The fundamental novel point here is the assumption of the thermal origin of the 
X and 7 radiation of the afterglow in the comoving frame of the shock front. The 
fit of the data in Figs. ^El gives a most clear and unambiguous support from the 
observations to this theoretical approach. 

All the works in the current literature tries to explain the afterglow emission by 
a very complex process implying magnetic fields, jet-like ejecta, e miss ion by a for- 
ward shock and a reverse shock (see e.g. Pirarii2l van Paradijs et al.l^, MeszaroJ^ 
and references therein). In our approach we evidence the existence of a much sim- 
pler process, directed forward, basically spherically symmetric and originating by 
a simple thermal emission in the comoving frame of the shock. 

W e are grateful to an anonymous referee for pointing out that Blinnikov et 
aim did argue that nonthermally looking GRB spectra can indeed be formed by 
a superposition of a set of thermal black body spectra with a temporal power- 
law evolution of the temperature. In our treatment not only time but also space 
integration on the EQTS takes place. This effect was explicitly omitted in the 
interesting paper of Blinnikov et al.l^ while their instantaneous GRB spectra are 
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Fig. 4. The time-integrated spectrum of the radiation observed in GRB 991216. The low energy 
part of the curve below 10 keV is fit by a power-law with index a = —1.05 and the high energy 
part above 500 kcV is fit by a power-law with an index /3 < —16. 



thermal, in our approach each instantaneous spectrum is derived from an infinite set 
of fohations of events on the EQTS, each one characterized by a different thermal 
spectrum in the comoving frame boosted by a different relativistic transformation 
obtained from the EOM. 

We emphasize that these results are extremely sensitive to the structure of the 
EQTS and to the theoretical assumptions adopted for each GRB era (see examples 
in Ruffini et al."^, Bianco & Ruffim^^^^). Due to the enormous redundancy built 
into the almost 10^ different paths mentioned above, possibly unprecedented in 
physics and astrophysics, we can assert the uniqueness of the solution. We also 
conclude that there is a marked difference (see Fig. ^ between the bolometric 
intensity of the afterglow, with a simple power-law behavior with an index n = — 1.6 
in the decreasing part, and the actual luminosity in a fixed bandwidth, which can 
have a complex dependence on time. Such a complex behavior could be erroneously 
interpreted as a broken power-law supporting the existence of jet-like structures in 
GRBs. 

The physical reasons justifying the assumptions in Eqs.QH^I are presented in 
Ruffini et al!^. 
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